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Abstract 

A measurement of the forward-backward asymmetries of e + e~ — > bb and e + e~ — > cc events using 
electrons and muons produced in semileptonic decays of bottom and charm hadrons is presented. 
The outputs of two neural networks designed to identify b — > t~ and c — > £ + decays are used in a 
maximum likelihood fit to a sample of events containing one or two identified leptons. The b and c 
quark forward- backward asymmetries at three centre-of-mass energies y/s and the average B mixing 
parameter y are determined simultaneously in the fit. Using all data collected by OPAL near the Z 
resonance, the asymmetries are measured to be: 

A*>|= ( 4.7 ±1.8 ±0.1 )% Af B = (-6.8 ±2.5 ±0.9 )% at (^i) = 89.51 GeV, 
A%= ( 9.72 ±0.42 ±0.15)% Af B = ( 5.68 ± 0.54 ± 0.39) % at {y/s} = 91.25 GeV, 
A]* = ( 10 3 ±L5 ±o.2)% Af B = ( 14.6 ±2.0 ±0.8)% at (y/s) = 92.95 GeV. 

For the average B mixing parameter, a value of: 

X= (13.12 ±0.49 ±0.42)% 

is obtained. In each case the first uncertainty is statistical and the second systematic. These results 
are combined with other OPAL measurements of the b and c forward-backward asymmetries, and used 
to derive a value for the effective electroweak mixing angle for leptons sin 2 of 0.23238 ± 0.00052. 
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1 Introduction 



The measurement of the forward-backward asymmetries of heavy quarks, A^ (q=b,c), in e + e~ — ► qq 
events provides an important test of the Standard Model. The bb forward-backward asymmetry 
provides one of the most precise determinations of the effective electroweak mixing angle for leptons 
sin 2 9^ s (assuming lepton universality). This is of particular interest at the present time in view of 
the nearly three standard deviation difference between the average values of sin 2 9^ s derived from 
quark forward-backward asymmetries at LEP on the one hand and from lepton forward-backward 
asymmetries at LEP and the left-right asymmetry at SLD on the other pQ. 

The forward-backward asymmetry arises from the cos 9 term in the differential cross-section for 
e + e~ — ► qq, 

^_ K l + C08 2 fl+ | A q| COS ; (1) 
d cos 9 3 

where 9 denotes the angle between the outgoing quark and the incoming electron flight directions, and 
where initial and final state radiation, quark mass effects and higher order terms have been neglected. 
The asymmetry Ap^ is related to the vector, gy, and axial- vector, g\, couplings of the Z to the 
electron e and quark q. At the peak of the Z resonance and for the s-channel Z exchange process only, 
the pole asymmetry is given by 

,0,q = 3 . . _ 3 2ff Ve /ffAc 2ffy q / ff A q 

- 4 - 4 1 + {gVe/gAe) 2 ! + ( W5Aq )2 • W 

In the Standard Model the couplings for any fermion f are related to the fermion charge Qf and its 
effective electroweak mixing angle sin 2 9l s as follows: 

^ = (l-4|Q f |sm 2 ^ ff ). (3) 
9Af 

The values of sin 2 9^ s are all close to 0.25, so the value of the asymmetry parameter for electrons, 
A e , is small, and varies rapidly with sin 2 6^, but the value of A^ is large, approximately 0.94, and 
varies only slowly with sin 2 9^ s . This results in a relatively large forward-backward asymmetry for bb 
events, which is then very sensitive to sin 2 9^ via A e . 

This analysis uses hadronic Z decays observed by the OPAL detector at LEP to measure A F ^ and 
^4p c B . The event thrust axis is used to estimate the primary quark direction. Leptons produced in 
semileptonic decays of b and c hadrons, usually referred to as prompt leptons, are used to identify 
heavy flavour events, and their charge is used to distinguish between decaying quarks and antiquarks. 
The asymmetry for bb events is diluted by the effect of neutral B meson mixing [5] . This is quantified 
by the average mixing parameter, x, which is the probability that a produced b hadron decays as its 
antiparticle. The effect of B mixing is to reduce the observed asymmetry of bb events by a factor 
(1 — 2x) with respect to the asymmetry without mixing. The parameter x can be measured from the 
fraction of like-sign lepton pairs in events with an identified lepton in each hemisphere of the event. 
The b and c asymmetries, and the mixing parameter x-> are therefore measured in a simultaneous fit to 
events with one or two identified leptons. The contributions of bb and cc events to the lepton samples 
are separated from each other and from background by using several kinematic variables describing 
the lepton and its associated jet, combined using neural network flavour separation algorithms. 

This paper supersedes our previous publication [3]. Compared to the previous paper, this analysis 
benefits from the inclusion of data recorded after 1994, and from a reprocessing of the full OPAL 
data set with final tracking algorithms and detector calibrations, which have improved in particular 
the performance of the electron identification and flavour separation algorithms. The systematic 
uncertainties have been reduced due to better knowledge of heavy flavour production and decay, and 
several details of the fit have also been improved. Similar analyses have been published by the other 
LEP experiments Analyses identifying b quark events via the resolvable decay length of b 

hadrons 0|BliniE3 also give precise measurements of the b quark forward-backward asymmetry. In 
addition, the charm asymmetry has been measured using reconstructed charm hadrons |111I121[T3] . 
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2 Data sample and event simulation 



The OPAL detector is described in detail elsewhere |14|H5| llf)]. The central tracking system is located 
inside a solenoid which provides a uniform axial magnetic field of 0.435 T along the beam axis. The 
beam axis coincides with the z-axis of the detector, with the polar angle 6 measured with respect to 
the direction of the electron beam. The azimuthal angle <j) is measured around the z-axis, and the 
radius r is the distance from the z-axis. The innermost part of the tracking system is a two-layer 
silicon micro-vertex detector. The silicon detector is surrounded by a vertex drift chamber, a large 
volume jet chamber with 159 layers of axial anode wires and a set of z chambers measuring the track 
coordinates along the beam direction. The jet chamber also provides a measurement of the rate of 
energy deposition along the track, dE/dx. 

Outside the magnet coil, the tracking system is surrounded by a lead-glass electromagnetic calorime- 
ter. The barrel region covers the polar angle range | cos#| < 0.82, and full acceptance to | cos6>| < 0.98 
is provided by the endcaps. The iron return yoke and pole pieces of the magnet are instrumented with 
streamer tubes and thin multiwire chambers to act as a hadronic calorimeter. The calorimeters are 
surrounded by muon chambers. 

Some flavour separation aspects of this analysis exploit the long lifetimes of b and c hadrons, and 
therefore rely in particular on the silicon microvertex detector. This detector was first operational in 
1991, providing measurements in the r-<p plane only. In 1993 it was upgraded to measure tracks in 
both r-(j) and r-z planes [E|> and in 1996 the cos 6 coverage for at least one silicon measurement was 
extended from | cos^j < 0.83 to | cos#| < 0.93 |16j . The changing detector performance with time was 
taken into account in the analysis. The resolution of the vertex drift chamber is sufficient to ensure 
some flavour separation even in data taken without silicon microvertex detector information. 

Hadronic Z decays were selected using standard criteria, as in ^7j- All the LEP1 data collected 
at centre-of-mass energies close to the peak of the Z resonance between 1990 and 1995 were included, 
together with Z-peak data recorded for detector calibration purposes during the higher energy LEP2 
running between 1996 and 2000. The thrust axis direction was calculated using charged particle tracks 
and electromagnetic calorimeter clusters not associated to any track. The polar angle of the thrust 
axis #t was required to satisfy | cos#t| < 0.95. Over four million hadronic events were selected. The 
exact numbers are listed in Tabled As indicated in the table, some events were recorded at centre- 
of-mass energies above and below the Z peak. The bulk of these off-peak data were recorded in 1993 
and 1995 at centre-of-mass energies approximately 1.8 GeV away from the peak. The other off-peak 
centre-of-mass energy samples from 1990 and 1991 are combined with these main samples, yielding 
measurements of the forward-backward asymmetries at three separate energy points. 



Mean Energy 
of Z events 


Selected 
Z decays 


Single lept 
Electrons 


on events 
Muons 


Dilepton 
events 


peak-2 (^/i) = 89.51 GeV 
peak = 91.25 GeV 
peak+2 (y^) = 92.95 GeV 


194 211 
4 079 047 
278 257 


11 567 
239 505 
16 977 


19 809 
410 877 
30 066 


321 
6 014 
394 



Table 1: The number of events selected below, on and above the Z peak. The average centre-of-mass 
energies of the selected Z events, the numbers of Z events, and the numbers of events with identified 
leptons are shown. Note that dilepton events are only selected in a region of high b-purity. 

Monte Carlo simulated events were generated using JETSET 7.4 ^Hj with parameters tuned by 
OPAL ^]. The fragmentation function of Peterson et al. |20j was used to describe the fragmentation 
of b and c quarks. The semileptonic decay model of Altarelli et al. [2j| with parameters fixed by 
CLEO, DEL CO and MARKIII data [2lfll23U24j was used to predict the lepton momentum in the rest 
frame of b and c hadrons. The generated events were passed through a program that simulated the 
response of the OPAL detector [25] and through the same reconstruction algorithms as the data. 
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3 Lepton identification and flavour separation 



Leptons were identified in hadronic events using well established algorithms JJj. All tracks with 
momentum greater than 2 GeV and | cos#| < 0.96 were considered as lepton candidates. 

Electron candidate tracks were required to have at least 20 measurements of dE/dx from jet 
chamber hits. Electrons were then identified using a neural network algorithm. The identification 
relies on ionisation energy loss (dE/dx) measured in the tracking chamber, together with spatial and 
energy-momentum (E/p) matching between tracking and calorimetry. The neural network output 
was in the range to 1, and was required to be greater than 0.9 for electron candidates. Photon 
conversions were rejected using another neural network algorithm ^7]. The efficiency of this selection 
for genuine electrons, not considering electrons from photon conversions, is 66%. The purity is 73%, 
where photon conversions are included in the background. 

Muons were identified by requiring a spatial match between a track reconstructed in the tracking 
detectors and a track segment reconstructed in the external muon chambers. Further rejection of 
kaons was achieved with loose dE/dx cuts |17j . These requirements selected muons with an efficiency 
of 74 % and a purity of 53 %. 

In this analysis, only leptons from the decay of a b or c hadron in a primary bb or cc event are 
considered as signal. Any other genuine electron or muon, and any hadron misidentified as a lepton, 
is considered as background. With this definition, about 11% of the electron background is from 
genuine electrons in light quark events, 49 % from misidentified hadrons in any primary quark flavour 
event, and 40% from photon conversions. About 68% of the muon background originates from pions, 
and 28 % from kaons. 

The relationship between the lepton charge and the primary quark or antiquark from whose decay 
it originated is vital for the asymmetry measurement. Leptons coming directly from the weak decay 
of b hadrons are denoted b — > t~ . A negatively-charged lepton comes from the decay of a hadron 
containing a b quark, and a positive lepton from a b antiquark 1 . Electrons and muons from leptonic r 
decays where the r lepton comes from a direct b decay, b —* t~ — > £~ , have the same sign correlation 
as the b — > £~ events. Both lepton charges are possible if a b hadron decays to a c hadron, which 
then decays semileptonically. These decays are written as b — > c — > £ + and b — > c — ► £~ , according to 
the charm content of the intermediate meson, and are both called cascade decays. In this analysis, 
any identified leptons from cc mesons, e.g. J/ip — * £ + £~ decays, are included with the cascades of 
the appropriate sign. A neutral B meson may have mixed before decay, so that a primary b quark 
decays as a b antiquark, or vice versa. Leptons from these mixed mesons are classified according to the 
decaying b quark, and contribute to the asymmetry with the wrong sign for their category, causing a 
reduction in the observed bb asymmetry by a factor of (1 — 2y). Leptons from the decay of c hadrons 
produced in cc events, c — ► £ + , have the opposite sign correlation to direct b — ► £~ decays; a primary 
c quark gives a positive lepton. 

As in the previous analysis two neural networks |26| denoted NETb and NETc were used to 
separate b — ► £~ and c — ► £ + decays from each other, from the cascade decays, b — > c — > £ + and 
b —* c — » £~ which dilute the observed forward-backward asymmetries, and from backgrounds. The 
networks were retrained to take full advantage of the improved detector calibrations. Several of the 
network input variables refer to the jet containing the lepton track. The same tracks and clusters used 
to define the event thrust axis were combined into jets using a cone algorithm |27| with a cone half- 
angle of 0.55 rad and a minimum jet energy of 5.0 GeV. The transverse momentum, pt, of each track 
was defined relative to the axis of the jet containing it, where the jet axis was calculated excluding 
the momentum of the track. A lepton sub-jet was defined as in jHJ. The sub-jet includes particles that 
are nearer to the lepton than to the jet axis, and is a measure of the lepton isolation. 

The first network, NETb, was trained to distinguish between b — > £~ events and all other categories. 
The input variables were the momentum p and transverse momentum pt of the lepton candidate, the 

1 Charge conjugate decays are implied throughout this paper. 
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energy of the lepton sub-jet, 2? su b-j e t, the total visible energy of the jet (calculated using all tracks 
and unassociated calorimeter clusters comprising the jet), EJ^, and the scalar sum of the transverse 
momentum of all tracks within the jet, (X^Pt)jet- Separate networks were trained for electrons and for 
muons. For the electron nets, two extra variables were included, namely the outputs of the electron 
and conversion neural networks. The distributions of the NETb output in the OPAL data are shown 
in Figure E^a) and (b), together with the predictions from Monte Carlo; at high NETb values, the 
separation of b — * £~ events from all other lepton sources is clearly visible. For example, requiring 
NETb > 0.7 gives a sample 89 % pure in b — > £~ decays, with 3 % b — ► c — > £ + and 4 % c — > £ + , whilst 
retaining 45% of all b — > t~ decays with an identified lepton (averaged over electrons and muons). 
The Monte Carlo gives a reasonable description of the data, and the effect of the small discrepancies 
visible is discussed in Section IB~31 

The second network, NETc, was trained to distinguish c — > £ + events from all other categories, 
including b — > t~ . The network used all the NETb variables, including the electron and conversion 
neutral network outputs for electron candidates, together with the following three quantities: the decay 
length significances, (L/ctl)i,2, of the jet containing the lepton (jet 1) and of the most energetic of the 
other jets in the event (jet 2), and the impact parameter significance of the lepton with respect to the 
primary vertex, d/a^. The decay length significance is the distance between the primary vertex and a 
secondary vertex, reconstructed from a subset of the tracks in the jet using the algorithms described 
in divided by the error on the decay length. The impact parameter significance is the distance 
of closest approach of a track to the primary vertex, divided by its error. The distributions of the 
NETc output in the OPAL data are shown in Figure W[c) and (d), together with the predictions from 
Monte Carlo; again the separation of c — ► £ + events from other leptons is clearly visible. Requiring 
NETc > 0.7 gives a sample 59 % pure in c — > £ + decays, with 12 % b — > £~and 6 % b — ► c — > £ + , whilst 
retaining 19 % of all c — » £ + decays with an identified lepton. 

The values of NETb and NETc were evaluated for all lepton candidates. When more than one 
lepton was identified in an event, the candidates were ranked according to the value of NETb. If the 
best two lepton candidates in an event satisfied NETb > 0.6, and were in opposite thrust hemispheres, 
then both candidates were retained, and this was classified as a dilepton event. Otherwise only the best 
lepton candidate was considered, and the event was classified as a single lepton event. The numbers 
of single and dilepton events selected are given in Tabled 



4 Fit method and results 

The values of the asymmetries and x were determined using a simultaneous fit to the observed numbers 
of single lepton events as a function of cos #x> NETb and NETc, and the numbers of dilepton events as 
a function of cos#t- The thrust axis of each event was used to estimate the direction of the primary 
quark, and the charge Qi of the identified leptons was used to distinguish between the quark and 
antiquark direction. 

The total likelihood to be maximized is the product: 

— ^single e X /2 s ingle fi X ^-double • (4) 

The single and double lepton likelihood terms are discussed in more detail below. 



4.1 Single lepton likelihood 

For single lepton events, the observable y = — Q^cos^t was used to classify events as forward, with 
y > 0, or backward, with y < 0. The observed forward-backward asymmetry was then examined in 
bins of NETb, NETc and \y\. The fit considers four classes of leptons according to the lepton charge 
and primary quark flavour: 

1. b -> £-, b -» r- -> £~ and b -» c -> tr 
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Figure 1: Outputs of the neural networks designed to select (a) b — ► e, (b) b — > fx, (c) c — > e and (d) 
c — ► /x decays. The data are shown by the points with error bars, and the expected contributions from 
different lepton sources by the hatched histograms. 
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4. background 



Direct b — ► t~ decays and the other decays in class (1) contribute to the observed asymmetry with 
the opposite sign to classes (2) and (3). The asymmetry of events coming from b decays is scaled by 
a factor (1 — 2%) for category (1) decays, and (1 — 2rjx) for category (2), to account for the effects of 
neutral B meson mixing. The additional correction factor T) takes into account that the two samples 
include different fractions of different species of b hadron: B°, B + , B^, b-baryons etc., due mainly 
to the different semileptonic branching ratios of D° and D + mesons. The value of rj is set to 1.083, 
evaluated from the Monte Carlo simulation. No equivalent correction factor was used for the small 
fraction of b — * c — > t~ events included in category (1). The fractions of each single lepton category 
in the full data sample are given in Table EJ Note that these are fitted fractions. The fit adjusts the 
overall rate of background in the sample. 





electrons 


muons 


(1) 


b- 


+ f 






35.6 % 


24.6 % 




b- 


-> c - 






2.8 % 


2.3 % 




b- 




-►r 




1.1 % 


0.7 % 


(2) 


b- 


-> c - 


+ i+ 




12.6 % 


9.7 % 


(3) 


c — 








19.3 % 


14.9 % 


(4) 


background 




28.6 % 


47.8 % 



Table 2: The composition of the single lepton samples. The first three lines are combined into category 
(1) in the fit. The fractions given correspond to the fractions j[ after the fit. 





ee 


MM 


e/i 


(1)(1) 


86.9 % 


86.3% 


86.4% 


(1)(2) 


10.3% 


8.5% 


9.6% 


(2)(2) 


0.3% 


0.2% 


0.3% 


(3)(3) 


0.4% 


0.4% 


0.5% 


(1)(4) 


1.8% 


4.0% 


2.8% 


others 


0.3% 


0.6% 


0.4% 


total number 


1 308 


2 067 


3 354 



Table 3: Fractions of events in the most important categories in the double lepton sample, together 
with the total number of ee, mi and efi events. 

The likelihood for single-lepton events has two terms: 



•^single I — Yl II 



(n F + n B )! /l+^FB\ np (I- A 



'FB 



"B 



NETb,NETc l-l " F< "° ! V 2 

x n -^expf-i^)^ (5) 

NETb.NETc VZ7rcTm V m / 

The first term is a product over bins of NETb, NETc and \y\ of binomial probabilities for the number 
of forward and backward events observed in the data in each bin, np and na, depending on the 
asymmetry Afb expected in the bin. The fit used 10 equally spaced bins in each of NETb, NETc and 



(Vi) (GeV) 



ss 

uu 

dd 



A FB ( %) (ZFITTER prediction) 
89.51 91.25 92.95 



5.95± 5.0 
-3.08±20.0 
5.96±10.0 



9.64±1.3 
6.32±7.3 
9.64±3.7 



11.89± 5.0 
12.07±20.0 
11.89±10.0 



Table 4: The values of the forward-backward asymmetries for light flavours taken from ZFITTER, 
and the variations used for calculating the background asymmetry. 
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muons 




j-qq 

J back 


qq 
c dilute 


j-qq 

J back 


qq 

C dilutc 


bb 


0.139 


0.028 


0.192 


0.093 


cc 


0.146 


-0.104 


0.178 


-0.043 


ss 


0.214 


0.025 


0.243 


0.103 


uu 


0.220 


-0.114 


0.169 


-0.123 


dd 


0.218 


0.107 


0.218 


0.083 



Table 5: The fractions of background coming from each primary quark flavour for electrons and 
muons, and the dilution factors by which the primary quark forward-backward asymmetry are scaled 
to evaluate the background asymmetry. 

\y\. The second term constrains the total number of events nx in a given NETb-NETc bin (with no 
binning in \y\) to the expected number, m, and allowed the overall normalisations of the electron and 
muon backgrounds to be determined from the data, reducing the overall uncertainty. 

The probability for an event to be forward is (1 + j4fb)/2 where A f -q is the expected forward- 
backward asymmetry in the bin considered: 

4 8 11 

A FB (NETb,NETc,|y|) = ]T /^(NETb, NETc, |y|),4 FB r:rfV (6) 

i=l ^ 

In this expression, f[ denotes the fraction of leptons in class i, and pi the normalised distribution of 
leptons from this class in bins of NETb, NETc and \y\, which is taken from Monte Carlo simulation. 
The source fractions j\ are derived from the Monte Carlo fractions /j. However, the fraction of 
background is a free parameter in the fit, and the fractions of prompt sources with i = 1,2,3 are 
given by f, t = (1 — fh)fi/(fi + fi + f"&)- In this way, the relative rates of the prompt leptons are fixed 
by the Monte Carlo simulation, and the fractions satisfy Ylifl = 1- The nominal asymmetries Ap B 
for each class are: 

A^ B = (l-2xK b B _ 

A 2 FB = - (l-2r?x)A^ (7) 

A3 . . _ acc 
^FB — ^FB 

.4 a Background 

A pB — Ji FB 

The background asymmetry J 4^ ck § round a j so depends on the primary quark asymmetries. It is ex- 
pressed as a sum over quark flavours: 

^.Background _ f<M r m A<M (o\ 
^FB — ^back c dilute /1 FB ■ \°J 

q 

This introduces a very weak additional dependence on the fitted values of the b and c asymmetries. 
The central values of the light quark asymmetries are taken from the predictions of ZFITTER 6.36 |28j . 
and are listed in Table |1J The quoted uncertainties in these asymmetries are taken from measurements 
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of the strange quark asymmetry by DELPHI |29j and of the light quark asymmetries by OPAL |30| . 
These measurements are consistent with the ZFITTER expectations, with relatively large statistical 
errors. The fractions /^L of each quark flavour contributing to the background are taken from 
the Monte Carlo prediction, as are the dilution factors Cj? lute which take into account the fraction 
of the primary quark asymmetry that is seen in background events of this flavour. The fractions 
and dilutions are listed in Table [5J Some contributions to the background, for example for photon 
conversions, have zero forward-backward asymmetry. Others, in particular kaons in ss events, inherit 
a significant fraction of the primary quark asymmetry. The background asymmetry actually varies as 
a function of NETb and NETc, but this effect is neglected in the fit as discussed in Section f5 . 31 below . 

In the second term in the single lepton likelihood, the number of events, m, expected in a bin of 
NETb and NETc is calculated by: 

4 

m(NETb, NETc) = iv£ ata £ /^(NETb, NETc) , (9) 

i=i 

where N^ ata denotes the total number of single lepton events in the data (I = e, p), and pi is the 
normalised distribution in bins of NETb and NETc for leptons in class i. The observed number of 
events in a bin, nx, is assumed to be Gaussian distributed with mean m and standard deviation o m 
computed as = m + (6^ c ) 2 , where <5^J is the uncertainty on the expected number m of events in 
this bin due to the limited Monte Carlo statistics. 



4.2 Double lepton likelihood 

The double lepton likelihood is a product in bins of \y\ of multinomial probabilities for the event 
to be forward, pp, backward, pb, or same-sign, p$. Both leptons have NETb > 0.6, and no further 
subdivision in NETb or NETc is made. The composition of the double lepton sample is given in 
Table 01 The ee, ep and pp events are all considered together: 

^double = I] (nF + nB | + 7 )! ( PF r (PET* (PST 5 ■ (10) 

, np\ rtR! ns- 

\y\ 

In this case, np, ne and ns are the number of forward, backward and same-sign events observed in 
the data in this \y\ bin. A dilepton event with opposite charge leptons is called forward or backward 
according to the direction of the negative lepton. The probability of such an event being forward or 
backward usually depends on the forward-backward asymmetry of the dilepton category, ij, where 
the two indices refer to the categories i = 1 — 4 of the two leptons. The fraction of same-sign dilepton 
events depends only on the mixing parameter x, and not on the asymmetries. Note that the fraction 
of forward and backward events with one lepton from category 1 and the other from category 2 is also 
independent of the forward-backward asymmetry, and only depends on the mixing parameter. This 
is because in the absence of mixing these events would all be same sign, losing the information about 
the asymmetry. If they are of opposite sign, without knowing which lepton comes from the mixed 
meson, the asymmetry information cannot be recovered. 

Only certain categories ij are possible. Any flavour event can give a background lepton, but lepton 
classes 1 and 2 can only come from bb events, and class 3 only from cc events. The possible dilepton 
classes are therefore: ij = 11, 12, 22, 33, 14, 24, 34, 44. The overall probabilities for forward, backward 
or same sign events are given by a sum over the possible classes ij: 

P x = J2M\y\)p%(\y\) (ii) 

Since the contribution of background leptons to the double-lepton sample is small (see Table EJ , the 
fractions are taken directly from the Monte Carlo simulation, without being corrected for the rates 
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of background leptons fitted in the single lepton sample. Similarly, the small residual asymmetries 
of the background leptons are neglected. Writing Y = 8|y|/3(l + y 2 ), the probabilities B s for an 
event to be forward, backward or same-sign for each dilepton category are given by: 
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4.3 Results 

The data are divided into three separate energy bins whose mean energies are shown in Tabled The 
asymmetries were fitted simultaneously at all energy points, the overall likelihood being the product of 
the likelihood in Equation0]for each point. Therefore, the fit has nine free parameters: the asymmetries 
Ap B and Ap c B at three energy points, together with values common to all three energy points for the 
mixing parameter x, and the electron and muon background fractions, /£. Data from all years were 
fitted simultaneously, with the various fractions and (NETb, NETc) distributions determined from an 
appropriate mix of Monte Carlo events with different simulated detector configurations, taking into 
account the changes in the performance of the OPAL detector over time. 

The result of the fit is illustrated in Figure |^1 For this plot, two regions of NETb-NETc were 
selected, one 93 % pure in b — > t~ events and the other 59 % pure in c — > £ + events. The asymmetry 
observed in the data in bins of \y\ and the predicted asymmetry according to Equation are shown. 
The predicted asymmetry is calculated with the fitted results. The sign of the observed asymmetry is 
clearly different in the two regions. 
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A small correction is applied to the fitted asymmetries to correct for the effects of gluon radiation 
from the primary quark pair and the approximation of the original quark direction by the experimen- 
tally measured thrust axis jHJ. The effects of gluon radiation have been calculated to second order 
in a s using the parton level thrust axis to define the asymmetry [22], and the translation from the 
parton level to the hadron level thrust axis (defined using all final state particles without detector 
effects) has been determined using Monte Carlo hadronisation models f2T|. However, the detector 
acceptance and fitting method introduce additional biases which reduce the sensitivity of the anal- 
ysis to these QCD corrections. Therefore, the final correction applied to the fitted asymmetries is 
determined by comparing the asymmetries fitted on large samples of Monte Carlo simulated events 
with the true quark level asymmetries, after scaling the hadron level thrust axis QCD corrections in 
the Monte Carlo to the theoretical values [SB1- This procedure also accounts for any other residual 
biases in the fit, for example from the treatment of the background asymmetry. Combining all ef- 
fects, the raw fitted asymmetries were scaled by factors of 1.0050 db 0.0063 ± 0.0050 for b quarks and 
1.0117 db 0.0063 db 0.0062 for c quarks to determine the quark level results, where the first errors result 
from theoretical uncertainties 31,33 and the second from limited Monte Carlo statistics. 

After correcting for QCD effects, the results are: 



Abb _ 

Abb _ 
^FB" 

Abb _ 
^FB" 



/l 



FB" 
acc _ 
^FB" 



4' 



( 4.7 ±1.8 ±0.1 )% 
( 9.72 ±0.42 ±0.15)% 
(10.3 ±1.5 ±0.2 )% 

For the average B mixing parameter, a value of: 

X= (13.12 ±0.49 ±0.42)% 



FB" 



(-6.8 ±2.5 ±0.9)% at (y/s) = 89.51 GeV, 
( 5.68 ±0.54 ±0.39)% at (y/s) = 91.25 GeV, 
( 14.6 ±2.0 ±0.8 )% at (s/s) = 92.95 GeV. 



is obtained. In each case the first error is statistical and the second systematic. The evaluation of 
the systematic errors is described in the following section. The statistical correlations between the 
results are given in Table El The background levels in the electron and muon samples were fitted to be 
(89.4 ±0.6) % and (94.4 ±0.3) % of the rates estimated from Monte Carlo simulation, where the quoted 
errors are statistical only. These values are consistent with the known uncertainties in modelling the 
lepton background levels [T7] . 
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0.17 


0.10 
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^FB 




1.00 


0.00 
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^FB 




1.00 


0.00 


X 






1.00 


X 






1.00 



Table 6: Statistical correlation matrices at each centre-of-mass energy using the entire (1990-2000) 
data sample. The correlations with the fitted background fractions are negligibly small (< 0.01) and 
are not given here. 



5 Systematic uncertainties 

Systematic uncertainties result from a number of sources, including modelling of b and c hadron 
production and decay, external branching ratio inputs, background uncertainties and the performance 
of the OPAL detector. For the first two sources, the proposals of the LEP Heavy Flavour Electroweak 
Working Group have been adopted [SSI- Where errors are only assessed by a variation in one direction, 
for example to an alternative model, or when slightly asymmetric errors are assessed, the larger 
deviation is taken as la. The systematic errors are summarised in Tables [7| and [S] and discussed in 
more detail below. The signs of the errors indicate the direction of change in the result when the 
corresponding quantity is varied. 
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Figure 2: An illustration of the fit results for (a) a b enriched region and (b) a c enriched region 
of NETb-NETc space. The asymmetry observed in the data in each bin of \y\ is compared with the 
expectation calculated from the full fit. The errors are purely statistical. 
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Af B (%) 


Af B (%) 


x(%) 


Fitted Value 


9.716 


5.683 


13.121 


Statistical error 


±0.418 


±0.542 


±0.485 


Systematic error 


±0.150 


±0.386 


±0.421 


Sources of systematic errors 


b — ► I semileptonic decay model 


±0.011 


T0.101 


±0.175 


c — > £ + semileptonic decay model 


±0.064 


T0.036 


T0.286 


(x h E ) ±0.008 


T0.012 


T0.008 


TO. 074 


\x%) ±0.008 


±0.054 


T0.024 


±0.010 


Total models 


0.085 


0.110 


0.344 


BR(b-»£ - ) (10.65 ±0.23)% 


=F0.020 


±0.104 


±0.120 


BR(b -> c ->■ £+) (8.08±0.18) % 


T0.006 


T0.050 


T0.105 


BR(b -► c £-) (1.62±0.44) % 


±0.011 


±0.174 


±0.014 


BR(b ->■ r- ->■ £-) (0.419±0.05) % 


to.ooi 


±0.028 


±0.007 


BR(c £+) (9.77±0.32) % 


±0.028 


T0.151 


±0.008 


# b (21.647±0.072) % 


T0.003 


±0.009 


0.000 


R c (16.83±0.47) % 


±0.019 


T0.099 


±0.002 


Total branching ratios 


0.041 


0.277 


0.161 


e ID efficiency ± 4.1 % 


=F0.005 


=F0.007 


±0.004 


H ID efficiency ± 3.0 % 


T0.002 


T0.021 


±0.022 


Conversions ±15 % 


±0.001 


±0.012 


T0.007 


Muon background (K) ±30 % 


±0.011 


±0.023 


TO. 034 


Muon background (ir) ±5 % 


±0.004 


±0.026 


T0.015 


Muon background (other) ±100% 


±0.009 


±0.017 


T0.058 


Background asymmetry u,d,s events 


±0.002 


±0.102 


±0.002 


Background flavour separation 


±0.020 


±0.020 


0.000 


Total background effects 


0.026 


0.114 


0.073 


Correction of % for b — > c — > £ + 


-0.008 


-0.069 


-0.068 


cos6*t dependence (fractions) 


-0.018 


±0.001 


±0.000 


Tracking resolution 


-0.025 


-0.046 


-0.016 


Flavour separation variables 


±0.067 


±0.163 


±0.099 


Time dependent mixing 


±0.032 


±0.098 


-0.091 


Monte Carlo statistics 


±0.050 


±0.041 


±0.070 


Charge reconstruction 


±0.015 


±0.004 


0.000 


LEP centre-of-mass energy 


±0.010 


±0.027 


0.000 


QCD correction 


±0.061 


±0.036 


0.000 


Total other systematics 


0.114 


0.216 


0.167 



Table 7: Results and breakdown of systematic uncertainties for the on-peak forward-backward asym- 
metries and the mixing parameter \- 
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±0.928 


±0.224 


±0.837 


Sources of systematic 


errors 






b — > £ semileptonic decay model 


_i_n oqi 
±U.Uol 


=F0.036 


=F0.020 


=F0.101 


c — ► semileptonic decay model 


-i-n mo 
^u.uiy 


±0.070 


±0.106 


T0.108 


\Xj^j ItU.UUo 


-|-U.UUo 


T0.037 


T0.019 


T0.034 


\Xfi) ±U.UUo 


-|-U.U1o 


±0.009 
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T0.012 
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±0.250 
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±0.003 


T0.002 


±0.003 


±0.025 


Muon background (K) ±30 % 


±0.003 


T0.032 


±0.010 


±0.073 


Muon background (tt) ±5 % 


±0.003 


T0.036 


±0.005 


±0.056 


Muon background (other) ±100% 


±0.002 


=F0.044 


±0.009 


±0.062 


Background asymmetry u,d,s events 


±0.018 


±0.832 


±0.006 


±0.302 


Background flavour separation 


±0.020 


±0.025 


±0.010 


±0.065 


Total background effects 


0.028 


0.836 


0.021 


0.333 


Correction of x f° r b ^ c ^ £ + 


-0.000 


-0.033 


-0.027 


-0.079 


cos #t dependence (fractions) 


±0.005 


-0.024 


-0.030 


±0.001 


Tracking resolution 


-0.011 


±0.261 


-0.002 


-0.215 


Flavour separation variables 


±0.029 


±0.172 


±0.097 


±0.403 


Time dependent mixing 


-0.037 


+ 0.036 


±0.069 


±0.112 


Monte Carlo statistics 


±0.058 


±0.102 


±0.064 


±0.118 


Charge reconstruction 


±0.007 


±0.005 


±0.016 


±0.010 


LEP centre-of-mass energy 


±0.007 


±0.018 


±0.003 


±0.007 


QCD correction 


±0.030 


±0.043 


±0.062 


±0.092 


Total other systematics 


0.082 


0.337 


0.155 


0.500 



Table 8: Results and breakdown of systematic uncertainties for the off-peak forward-backward asym- 
metries. 
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5.1 Modelling of b and c production and decay 

Semileptonic decay models: A correct description of the lepton momentum spectra in the rest 
frame of decaying b and c hadrons is crucial for the flavour separation. The semileptonic decays 
of heavy hadrons were described by the free-quark model of Altarelli et al. (ACCMM) with 
its two free parameters fixed by fits to CLEO data [22] f° r b — * £~ decays and the combined 
measurements of DELCO |23] and MARK III |24j for c — > £ + decays [33]. For cascade decays 
b — > c — > £ + and b — > c — > £~ , the D momentum spectrum measured by CLEO [33] is combined 
with the c — > £ + model to generate the lepton momentum distribution. Uncertainties in the 
CLEO D momentum spectrum are negligible compared with the uncertainties in the b — * l~ 
and c — > £ + models. 

For b — > £~ decays, the systematic uncertainties were assessed by reweighting the rest frame 
momentum spectrum according to the form-factor model of Isgur et al. [31], which has no free 
parameters (ISGW). This spectrum is harder than the central ACCMM model. As an alternative, 
the same model was used with the fraction of D** mesons produced in b decays increased from 
11 % to 32 % to describe better CLEO data, denoted ISGW**. This gives a softer spectrum than 
the default model. Although the weights to go from model to model were evaluated using only 
B° and B + meson decays, these weights were then applied to all b-hadron decays. Variations in 
the c — ► £ + decay model were assessed by varying the free parameters of the ACCMM model to 
give harder or softer decay models, constrained by the DELCO and MARK III data. For both 
b — * £~ and c — > £ + decays, the sign of the error assigned indicates the variation observed with 
the harder alternative spectrum. 

Heavy quark fragmentation: The lepton momentum spectra depend on the energy distributions 
of the heavy hadrons produced in bb and cc events. The Monte Carlo events were reweighted 
so as to vary the average scaled energy (xe) of weakly decaying b hadrons in the range (xe) = 
0.702 ± 0.008 and of c hadrons in the range (x E ) = 0.484 ± 0.008 [HS]- The fragmentation 
functions of Peterson et al, Collins and Spiller, Kartvelishvili et al. and the Lund group |20U36| 
were each used as models to determine the event weights, and the largest observed variations 
were assigned as the systematic errors. For b fragmentation, the largest shifts were observed with 
the function of Collins and Spiller, and for c fragmentation with the function of Peterson et al. 

5.2 Branching ratios and partial widths 

Semileptonic branching ratios: The values and uncertainties used for the branching ratios of 
semileptonic b and c hadron decays are listed in Tables [7] and [S] jU[33] • 

Partial widths of the Z: The fractions of hadronic Z decays to bb and cc were varied according to 
the uncertainties in the LEP average values 0: R h = 0.21647±0.00072 and R c = 0.1683±0.0047. 
In each case, the fraction of Z hadronic decays to light quarks was adjusted to compensate the 
variation. 

5.3 Background uncertainties 

Lepton identification: The analysis is only weakly sensitive to the lepton identification efficiency, 
since this mainly affects the ratio of prompt to background leptons, which is effectively deter- 
mined from the data via the fit of the background level. However, a small component of the 
background is composed of genuine leptons, so the lepton identification efficiencies were varied 
by ±4.1% for electrons and ±3.0% for muons 17 . 

Background composition: The overall background fraction is a fitted parameter. However, the 
background shape as a function of NETb, NETc and \y\ may be different for different contribu- 
tions to the background. The rate of untagged photon conversions was varied by ±15% |17| . In 
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addition, the \y\ distribution of conversions was varied according to the difference in the rate of 
identified conversions observed between data and Monte Carlo simulation. The contributions to 
the muon background from it, K and other sources were varied in turn by 5%, 30% and 100%. 
These uncertainties were evaluated using control samples of Kg — ► 7r + 7r~, three-prong tau decays 
and tracks with an enhanced kaon fraction selected using dE jdx cuts |37| . The variation of the 
source fractions, including backgrounds, as a function of cos 6<r is discussed in Section T5 .41 below . 

Background asymmetries: The uncertainties in the background asymmetry were evaluated by 
varying the light primary quark forward-backward asymmetries by the uncertainties listed in 
Table |U Note that the up and down asymmetries measured by OPAL are +91 % correlated 30 . 
They were therefore varied at the same time, and the resulting uncertainty combined in quadra- 
ture with the uncertainty in the strange asymmetry. The bb and cc asymmetries are treated 
self-consistently in the fit, and the uncertainty in these quantities is therefore automatically 
reflected in a very small contribution to the statistical uncertainty. 

Background flavour separation: The dependence of the background asymmetry on NETb and 
NETc is neglected in the fit (see Equation [8J) but any residual bias is removed by the bias 
correction as discussed in Section f4.3[ providing the Monte Carlo gives an adequate description 
of the variations. To quantify the size of any mismodelling, an additional fit was performed using 
a binned background asymmetry without changing the bias correction, and half the difference 
between this and the standard fit result was assigned as a systematic error. 

5.4 Other uncertainties 

Correction factor for mixing in cascade decays: The difference between using the standard value 
of rj = 1.083 and setting it equal to 1.000 was assigned as a conservative estimate of the uncer- 
tainty. 

Dependence of source fractions on cos 6^: The description of the source fractions as a function 
of cos #t was checked by selecting b — > £~ , c — > i + and background enriched regions of the NETb- 
NETc plane, and comparing the cos#t distributions of data and Monte Carlo simulation. No 
systematic trend was seen for the c — * £ + enriched sample. The ratio for background varied 
by up to 10 % in the endcap region, with smaller deviations for the b — > l~ enriched sample. 
A systematic uncertainty was therefore assessed by reweighting the Monte Carlo background 
according to the data/Monte Carlo ratio observed for the background enriched sample, and 
repeating the fit. The change in the background sample was compensated by an increase in the 
prompt leptons, so this procedure also accounted for any possible discrepancy in the b —* l~ 
sample. 

Tracking resolution: A global 10 % degradation in the tracking resolution was applied to the Monte 
Carlo sample, separately in the r-(j> and r-z planes, as discussed in [T7HTU] . The fit was repeated 
with each of these modified samples, and the sum in quadrature of the two shifts in the results 
was taken as a systematic uncertainty. 

Flavour separation variables: Possible uncertainties due to mismodelling of the neural network 
input variables were taken into account by reweighting the Monte Carlo events by the ratio 
between the data and Monte Carlo distributions of each variable in turn, and repeating the fit. 
The sum in quadrature of the shifts in the results was assigned as a systematic uncertainty, 
dominated by the result of reweighting the lepton p and pt distributions for the asymmetries, 
and the input parameter significance for x- The neural network output distributions were also 
reweighted, to account for the discrepancies between data and Monte Carlo visible in Figure ^ 
but the resulting changes in the results were much smaller than those seen when reweighting the 
inputs, and no additional systematic error was assigned. 
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Time dependent mixing: The mixing parameter x reflects the fraction of mixed events in an un- 
biased sample of all b hadrons. However, the time dependent oscillations for B° mesons are 
sufficiently slow that the lifetime variables used in NETc might change the effective value of x 
as a function of NETc. To evaluate the impact of such an effect, the value of x was evaluated 
for all bb Monte Carlo events with a lepton candidate, in bins of NETc, and for the entire bb 
sample. The fitted value of x in the likelihood function was then scaled by the ratio of Monte 
Carlo values of x m the NETc bin and x m an bb events. The shift in the fit results with this 
rescaling was taken systematic uncertainty. 

Monte Carlo statistics: The fit to the data was repeated 1000 times, with the Monte Carlo reference 
distributions being randomly varied according to the statistical uncertainty in each bin of NETb, 
NETc and \y\. The RMS variation in the results was added in quadrature with the Monte Carlo 
statistical uncertainty on the bias corrections discussed in Section fOl to give the total uncertainty 
due to limited Monte Carlo statistics. 

Lepton charge reconstruction and asymmetry: Monte Carlo simulation predicts that about 0.03 % 
of electrons and 0.3 % of muons are reconstructed in the tracking chamber with the wrong charge 
(electron tracks are subject to tighter quality requirements because of the stronger dE/dx cuts). 
This effect is accounted for in the fit, and the full size of the correction is taken as a systematic 
error. Studies in [SB] show a possible difference of around 10 -3 between the lepton identification 
efficiencies for positive and negative tracks. When combined with the small difference (2 %) in 
numbers of reconstructed leptons in positive and negative hemispheres, this leads to negligible 
asymmetry biases of around 2 x 10 -5 . 

LEP centre-of-mass energy: The LEP centre of mass energy is known to a precision varying be- 
tween 3.5 and 20 MeV depending on the energy point and year [22] • Taking year-to-year correla- 
tions into account, and assuming the Standard Model dependencies of the asymmetries on y/s, 
this leads to the uncertainties given in Tables [7] and [H] on the asymmetries at the quoted values 

of y^S. 

QCD and thrust axis correction: The corrections applied to the raw fitted asymmetries are known 
to precisions of 0.80% for b quarks and 0.89 % for c quarks, as discussed in Section^ The theo- 
retical part of this uncertainty is assigned as the error on the QCD correction, whilst the Monte 
Carlo statistical uncertainty is accounted for separately as discussed above. 

5.5 Consistency checks 

The analysis was performed separately for data from 1990 and 1991, for each year separately from 
1992 to 1995, and for data from 1996 to 2000. The results were consistent among themselves and with 
the standard result from fitting all data simultaneously. Taking into account statistical errors only, 
the x 2 /dof value for the six measurements to be consistent with the same value was 5.9/5. The 
equivalent x 2 /dof values for A^ B and x were 2.4/5 and 5.3/5. Separate fits for electron and muon 
events were also made and showed similar consistency. The number of bins in y, NETb and NETc were 
varied, and again consistent results were obtained. To verify the correctness of the fit method itself, it 
was tested on large samples of test Monte Carlo samples generated with various b and c asymmetries 
and with source fractions and input NETb and NETc distributions generated to correspond to those 
in the full simulation. In these tests, the fit was found to be unbiased and to return correct statistical 
errors on the fitted parameters. 
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6 Conclusions 



The forward-backward asymmetries of e + e~ — ► bb and e + e~ — > cc events have been measured at 
three energy points around the Z peak, using electrons and muons produced in semileptonic decays of 
bottom and charm hadrons. The results, corrected to the primary quark level, are: 

Af B = ( 4.7 ±1.8 ±0.1 )% Af B = (-6.8 ±2.5 ±0.9 )% at (s/s) = 89.51 GeV, 
A%= ( 9.72 ±0.42 ±0.15)% Af B = ( 5.68 ± 0.54 ± 0.39) % at (a/s) = 91.25 GeV, 
Af B = (10.3 ±1.5 ±0.2)% Af B = ( 14.6 ±2.0 ±0.8)% at (y/s) = 92.95 GeV. 

For the average B mixing parameter, a value of: 

X= (13.12 ±0.49 ±0.42)% 

is obtained. In each case the first error is statistical and the second systematic. The asymmetry results 
are shown clS EL function of a/5 in Figure 01 

Using the ZFITTER prediction for the dependence of A^ and A^ B on a/5, the measurements 
are shifted to mz (91.19 GeV), averaged and corrected for initial state radiation, 7 exchange, 7 — Z 
interference and quark mass effects. The results for the pole asymmetries are: 

A^ B = (9.81 ±0.40 ±0.15)%, 
A% B = (6.29 ± 0.52 ± 0.38) % . 

These results are consistent with those of A^ derived from inclusive jet charge measurements JO] 
and those of and A^ B derived from fully reconstructed charm hadrons ^H] • Whilst the statistical 
correlations between the asymmetries measured here and those from reconstructed charm hadrons are 
negligible, the jet charges of some of the events involving b —* l~ decays are also used to measure A^ 
in ^Hj. This results in a statistical correlation of 0.11 ±0.12 between the two measurements, evaluated 
using large samples of Monte Carlo simulated events. Taking both statistical and systematic correla- 
tions into account, and making very small corrections to the inclusive and charm meson asymmetries 
to use the same nominal centre-of-mass energies as in this analysis, the final OPAL results for the b 
and c quark asymmetries around the Z resonance are: 

A%= ( 5.3 ±1.2 ±0.1 )% Af B = (-4.8 ±2.2 ±0.7 )% at (^5) = 89.51 GeV, 
A%= ( 9.69 ±0.29 ±0.13)% Af B = ( 5.83 ± 0.49 ± 0.32) % at (a/5) = 91.25 GeV, 
A%= (11.4 ±1.0 ±0.2 )% Af B = ( 14.7 ± 1.8 ± 0.7 ) % at (a/s) = 92.95 GeV, 

and the pole asymmetries are measured to be: 

A°^ B = (9.89 ±0.27 ±0.13)%, 
A^ B = (6.57 ±0.47 ±0.32)%. 

In all cases, the first error is statistical and the second systematic. The x 2 /dof of the fit to all OPAL 
asymmetry measurements to give the two pole asymmetries is 11.2/13. (There are six asymmetry 
measurements from the lepton tag analysis, three from the inclusive analysis, and six from the charm 
meson analysis, which consists of a simultaneous fit to the b and c quark asymmetries.) Taking 
into account statistical and systematic uncertainties, these average pole asymmetries have a +15% 
correlation. Within the framework of the Standard Model, the measurements of Ap B and Ap B taken 
together can be interpreted as a measurement of the effective weak mixing angle for leptons of 

sin 2 6 e cS = 0.23238 ± 0.00052 . 

As can be seen from Figure |3 this result favours large values of the Higgs mass, in agreement with 
other measurements of heavy flavour asymmetries and in contrast to measurements of the leptonic 
forward-backward and left-right asymmetries 
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Figure 3: The measured b and c quark asymmetries as a function of centre-of-mass energy, for the 
lepton, inclusive (b only) ^U] and D meson (c only) analyses. The measurement of the b quark 
asymmetry from D mesons is of low precision and is not shown. The measurements from different 
analyses are slightly displaced on the horizontal axis for clarity. The combination of the lepton analysis 
and the other analyses is also shown, together with the Standard Model expectation for Higgs masses 
between 114 and 1000 GeV calculated using ZFITTER [H]. 
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